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Random matrix theory (RMT) has been used extensively and successfully in 
the field of nuclear physics, 1,2 where the predictions of the theory concerning the 
distribution of excitation energies and the distribution of transition strengths agree 
well with experimental data. The results of RMT concerning the excitation-energy 
distribution have been successfully applied also in condensed matter physics. 3-6 
In this work we demonstrate how RMT of transition strengths can be applied in 
condensed matter theory, in particular in the study of transport in the variable- 
range-hopping regime. We obtain universal behavior of the magnetoconductance 
(MC) in the deeply localized regime, which is determined by the crossover function 
between the orthogonal to the unitary ensemble, in the absence of spin-orbit scat- 
tering, and between the symplectic to the unitary ensemble in the presence of strong 
spin-orbit scattering. We derive this crossover function, and predict quantitatively 
the universal MC curves, which are verified numerically (Fig. 2). 

The idea underlying the application of RMT rests on the assumption that the 
statistical behavior of a complicated system is determined by its symmetries. Ac- 
cordingly, the Hamiltonian describing the system can be classified into one of three 
universality classes, characterized by a parameter /3 which counts the number of 
degrees of freedom associated with each matrix element: The Gaussian orthogonal 
ensemble {j3 = 1), the Gaussian unitary ensemble {(5 = 2), where time- reversal 
symmetry is broken, e.g. by a magnetic field, and the Gaussian symplectic ensem- 
ble (/3 = 4), where rotational symmetry is broken, e.g. by spin-orbit interactions. 
For these different ensembles one can derive the distributions of level spacings and 
of overlap probabilities, which in the case of nuclear transitions are related to the 
distributions of excitation energies and of transition strengths, respectively. 

The use of RMT in condensed matter problems has been confined so far to 
predictions concerning the level distribution. It was used for studying electronic 
properties of small metallic particles, 3 and it has been also invoked 4,5 to explain the 
theoretically predicted 7 universal conductance fluctuations in the weakly localized 
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regime in terms of the rigidity of the level spectrum of the transfer matrix 4 or the 
Hamiltonian. 5 Altshuler and Shklovskii 5 showed that when the Zeeman splitting is 
neglected, the amplitude of the conductance fluctuations is determined by a param- 
eter defined by x — 4(3 /s 2 , where s is the degeneracy of each level. Thus x — 
(a) 1, (b) 2, (c) 4, (d) 8 for the cases (a) both time-reversal and rotational symme- 
tries are conserved, (b) only time-reversal symmetry is broken, (c) only rotational 
symmetry is broken and (d) both symmetries are broken (f3 = 2). In particular, 
a magnetic field which breaks time-reversal symmetry decreases the magnitude of 
the fluctuations, independently of the amount of spin-orbit scattering present. This 
result was confirmed by experiments. 8 

More recently Pichard et al. 9 claimed that similar arguments can also be ap- 
plied to the strongly localized regime. Here the conductance is determined by an 
equivalent resistor network 10 in which each two impurities are connected by a con- 
ductance goJ(H) exp [— rij/£(H) — Ae/kT], where go has units of conductance, r\j 
is the distance between the impurities, Ae = (|ej| + |e 3 -| + |ej — e,-|)/2, where q are the 
energies of the impurity states, and £ is the localization length. Both £ and the am- 
plitude J depend upon the magnetic field. Using RMT predictions for the eigenvalue 
statistics, Pichard et al. concluded that £ is doubled as time-reversal symmetry is 
broken and increases by a factor of four as rotational symmetry is broken, in agree- 
ment with earlier exact results 11 in quasi-lci samples. 12,13 However, a relatively large 
magnetic field (= H^) of a unit quantum flux through an area £ 2 is necessary to 
induce a change in the localization length. (In such magnetic fields other effects, 
such as the shrinking of the wave functions may be significant, especially in doped 
semiconductors, where £ is on the order of the Bohr radius of the impurity state.) 
On the other hand the amplitude J(H) of the relevant hops is determined by the 
interference of all paths within a cigar shaped area of length R and width y^R^. 14,15 
Accordingly, the relevant magnetic field scale for a change in the amplitude, Hr, is 
a unit flux through the much larger cigar shaped area. (Experimental values for the 
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ratio between these two areas, (i?/£) 3 / 2 , 



range from 5 to 100 and more, depending on 



temperature 16 ). Consequently, the MC is determined by the amplitude, for a wide 
range of magnetic fields. Moreover, in the presence of strong spin-orbit scattering the 
localization length is unaffected by a magnetic field, 11 ' 13 and the MC is dominated 
by the magnetic field dependence of the amplitude. This amplitude is determined 
by the overlap between the impurity wavef unctions. The overlap distribution has 
been calculated using various numerical and analytic approximations. 14 ' 15 ' 17 ' 18 

In this work we use the fact that the overlap probability between two wave- 
functions can be written in a form analogous to the transition strength in nuclear 
physics, to demonstrate how the overlap distribution can be obtained similarly to 
the calculation of the transition-strength distribution, using RMT. Thus one expects 
this distribution to be of a universal nature, leading to universal predictions for the 
MC in the variable-range-hopping regime. 

In order to derive the crossover distribution we study a general interpolating 
N x N quaternionic matrix M, 



where Si(A ; ) are symmetric (antisymmetric) N x N random matrices, whose ele- 
ments are normally distributed with a zero mean and unit variance, and o\ are the 
Pauli matrices. For 7 = 5 = the matrix M belongs to the Gaussian orthogonal 
ensemble, for 7 = 1, 8 — it belongs to the Gaussian unitary ensemble, and for 
7 = 0, 8 = 1 it belongs to the Gaussian symplectic ensemble. Thus turning on 
the parameter 7 from zero to 1 interpolates between the Gaussian orthogonal and 
the Gaussian unitary ensembles when 8 = and between the Gaussian symplectic 
and the Gaussian unitary ensembles when 8—1. Similarly changing 8 from zero 
to 1, with 7 = 0, interpolates between the Gaussian orthogonal and the Gaussian 
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symplectic ensembles. M is normalized such that the average overlap is independent 
of the ensemble (i.e. of 7 and S). The overlap probability between sites i and j, 
summed over all final spin states and averaged over all initial spin states is given by 

y= l*rMtM^ = (l-£)(l-¥)^ + 
3 3 

J,^r^ o qx2 N „,2 r, r2 „,2 



i-^fE«I + (i-¥)^ + SE4 (2) 
fc=i fc=i 



where My is the 2x2 block of the matrix M, = (Sk)r,- and = (Ai)y. The 
distribution of the overlap probability, P(y), is readily expressed in terms of its 
Laplace transform F(s), 

POO 

F(s) = / e- sy P(y) 
Jo 

= V{[1 + 2(1 - - + 7 2 (1 - + f (1 - + ^] f } • 

For the four pure symmetry cases discussed above, which correspond to (a) 7 = 5 = 
0, (b) 7 = 1, 5 = 0, (c) 7 = 0, 5 = 1, and (d) 7 = 5 = 1, the overlap distribution 
function is given by 

p(v) = T^y x,2 ~ l ^ xy ' 2 - ( 4 ) 

Eq. (|4]) is the exact result for the transition-strength distribution for large matrices 
belonging to the Gaussian orthogonal (x = (3 = 1), the Gaussian unitary (x = 
(3 = 2) and the Gaussian symplectic (x = /3 = 4) ensembles. However, when both 
time-reversal and rotational symmetries are broken, P(y) is given by Eq. (|J), with 
X = 8, even though the 2N x 2N matrix M [Eq. (|I|)] belongs to the Gaussian 
unitary ensemble. The reason, as in the weakly localized regime, 5 is that even when 
time-reversal symmetry is broken, the broken rotational symmetry couples the two 
spin directions, so the overlap amplitude involves a sum over the two spins, leading 
to eight free variables [see Eq. (H)]. 

The distribution function can be expressed analytically for the crossover from 
the Gaussian orthogonal ensemble to the Gaussian unitary ensemble (5 = 0, finite 



7 ), p( y ) = e~yl z l^^fY^Z y jZ)l4Z , with Z = 2 7 2 (1 - 7 2 /2), which is the exact 
result for the transition-strength crossover function for this case, in the limit of large 
matrices. 2 ' 15 This function has been verified numerically for a kicked rotor. 19 Simi- 
larly, the crossover function from the Gaussian orthogonal ensemble to the Gaussian 
symplectic ensemble is given by 

1-1/2 



P(y) 



(i - ^-)*74 



y exp 



y(i-f) 

5 2 (l-f ) 



r 2/(1 -5 2 ) " 


-h 


\y(i-s 2 ) ' 




[<5 2 (l-^)_ 






i 



(5) 



while the breaking of time-reversal symmetry in the Gaussian symplectic ensemble 
(5=1) gives rise to the crossover function 
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For general 7 and/or 5, -P(y) can be written as a convolution of two functions. 
The small-?/ behavior, however, can be directly deduced from the Laplace transform 
F(s) and will be described by the power law corresponding to the lower symmetry 
ensemble, crossing over to the higher symmetry ensemble(s) behavior(s) at values 
of y which depend upon the values of 7 and S. 

The most important consequence of Eq. (f|) is that the power law that describes 
the small-?/ behavior, the region that contributes the most in the strongly localized 
regime, increases with a magnetic field (i.e. with 7), independently of the amount 
of spin-orbit scattering (i.e. of 5). This leads to a positive MC in the presence 
and in the absence of spin-orbit scattering. 13-15,17,18 We checked this prediction 
by calculating the distribution of the spin-averaged transmission probability T, at 
energy E/V = 0.1, through a 5x5 diamond, described by the Anderson Hamiltonian 
of band width 4V, with on-site uniform disorder of width W. We include an orbital 
magnetic field, characterized by the overall flux through the diamond 0, in units of 
0o, the quantum flux, and spin-orbit scattering, characterized by the typical angle 
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of rotation per hop in spin-space, A. 20 In Fig. 1 we plot the resulting small-T 
distributions for four cases (a) (p = 0, A = 0, (b) (p = 5.5, A = 0, (c) (p = 0, A = 2n 
and (d) (p = 5.5, A = 2n, corresponding to the four cases discussed above. The 
two panels in Fig. 1 correspond to the choices (a) W/V = 2 and (b) W/V = 4. 
We also plot in each panel the four slopes resulting from Eq. (d). There is clearly 
a satisfactory agreement between the predictions and the numerical results. In 
particular, there is a clear increase in the slope upon application of a magnetic field. 
As the localization length only determines the average transmission probability, the 
small-?/ behavior is insensitive to the value of disorder. 

In order to calculate the conductance one has to solve the percolation problem of 
the random resistor network, 21 ' 15 with resistors distributed according to P{y). This 
procedure leads to positive or negative MC for a system close to the metal-insulator 
transition, depending on the system parameters. 15 However, away from the metal- 
insulator transition, the conductance of the system is given 14 ' 18 ' 22 by G = Goe <ln<yyS}> , 
except for an exponentially small region of small magnetic fields. Go denotes a 
typical conductance of a single hop, while < ... > denotes an average over the 



distribution P(y). Thus the relative MC, 5G = [G(H) - G(H = 0)] /G(H = 0), will 



be universal — depending only on the change in P(y) as time-reversal symmetry 
is broken, namely as 7 increases from zero. In particular, this crossover function is 
given for the case of no spin-orbit scattering (5 = 0) by 5G = 71/2 — 7 2 . Since we 
expect the time-reversal breaking parameter 7 to be proportional to magnetic field, 
the MC is linear at small fields, 14,17,18 saturating at a value of 1 for a strong enough 
magnetic field. 15 For strong spin-orbit scattering (5 = 1) the relative MC is given by 



At low fields the MC is now quadratic, saturating at a value of exp(5/6)/2 — 1 ~ 0.15 
at high fields. 

We checked numerically these universal behaviors by calculating the MC through 
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the same 5x5 diamond giving rise to Fig. 1. The conductance for a single realization 
was calculated from the transmission probability, using Landauer formula, 23 and was 
logarithmically averaged over 100000 or 200000 realizations for W/V = 25, without 
spin-orbit scattering (Fig. 2a) and with strong spin-orbit scattering (Fig. 2b). The 
high disorder was chosen so that the localization length will be smaller than the size 
of the sample, giving rise to the separation of the magnetic field scale at which the 
localization length changes (H^) from the magnetic field scale at which we expect 
our predictions to hold (H R ). As is clearly seen in Fig. 2, there is a range of 
magnetic fields, corresponding approximately to one unit flux through the sample, 
where the MC, as expected, saturates. At higher fields, the MC increases further as 
the localization length starts to increase significantly. This effect is less important 
in the presence of spin-orbit scattering because here, for large enough samples, the 
localization length remains unchanged by the magnetic field. We also plot in the 
figure the prediction of our theory, in excellent agreement with the numerical data. 
The only fitting parameter in both panels was the ratio between the time-reversal 
symmetry breaking parameter, 7, and the magnetic flux, 0, which was chosen as 
7 = 1.50. The shape of the curves is very similar to what one expects in experiments: 
the MC at small fields should follow our universal predictions, while at higher fields 
other effects, such as the change in the localization length and the shrinking of 
the wavefunction start to play a significant role, and one expects deviations from 
universality. The deeper in the localized regime, the more separated the magnetic 
scales H R and become and the more universal the MC should be. 

To conclude we have derived the crossover distribution for the overlap probabil- 
ities (the transition strengths) between the different ensembles. This distribution 
reduces to all the exact random-matrix-theory results in the appropriate limits. This 
function has been used to predict the universal magnetoconductance curves deep in 
the localized regime, where the localization length is much smaller than the hopping 
length, and the percolation criterion for the conductance is equivalent to the log av- 
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eraging procedure. In that case there is a wide range of magnetic fields, where our 
universal predictions hold. In this regime we expect all the relative magnetoconduc- 
tance curves, e.g. for different temperatures, to collapse onto a single curve, when 
expressed in terms of the scaled time-reversal breaking parameter, presumably the 
magnetic flux through the cigar shaped hopping area. So far experiments have been 
confined to the vicinity of the metal-insulator transition, since in this regime the 
conductivity is more easily measurable. While one has to use more sensitive tools 
to probe the deep localized regime, we hope that our work will stimulate further 
experimental effort in this direction. 
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administered by the Israel Academy of Sciences and Humanities and the German- 
Israel Foundation for scientific research and development. 
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Figure Captions: 

1. The distribution of the transmission probability through a 5 x 5 diamond, de- 
scribed by an on-site disordered Anderson Hamiltonian, in the presence of magnetic 
flux and spin-orbit scattering for two values of disorder. The four curves correspond 
to (a) no magnetic field and no spin-orbit scattering, (b) finite magnetic field (c) 
finite spin-orbit scattering and (d) finite magnetic field and finite spin-orbit scatter- 
ing. Also depicted are the slopes (with arbitrary offsets) expected from Eq. (f|), (a) 
-1/2, (b) 0, (c) 1 and (d) 3. 

2. The magnetoconductance through a 5 x 5 diamond for disorder W/V = 25. 
The solid lines are the predictions of the theory, where the ratio 7/0 = 1.5, between 
the time-reversal symmetry breaking parameter 7 and the magnetic flux 0, is the 
only fitting parameter. 
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